Generation of magnetic field by dynamo action in a turbulent flow of liquid sodium 
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We report the observation of dynamo action in the VKS experiment, i.e., the generation of 
magnetic field by a strongly turbulent swirling flow of liquid sodium. Both mean and fluctuating 
parts of the field are studied. The dynamo threshold corresponds to a magnetic Reynolds number 
R m ~ 30. A mean magnetic field of order 40 G is observed 30 % above threshold at the flow 
lateral boundary. The rms fluctuations are larger than the corresponding mean value for two of the 
components. The scaling of the mean square magnetic field is compared to a prediction previously 
made for high Reynolds number flows. 
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The generation of electricity from mechanical work has 
been one of the main achievements of physics by the end 
of the XIX th century. In 1919, Larmor proposed that a 
similar process can generate the magnetic field of the sun 
from the motion of an electrically conducting fluid. How- 
ever, fluid dynamos are more complex than industrial 
ones and it is not easy to find laminar flow configurations 
that generate magnetic fields Two simple but clever 
examples have been found in the seventies J2J and have 
led more recently to successful experiments [3| . These ex- 
periments have shown that the observed thresholds are 
in good agreement with theoretical predictions Q made 
by considering only the mean flow, whereas the satura- 
tion level of the magnetic field cannot be described with 
a laminar flow model (without using an ad-hoc turbu- 
lent viscosity) [5j. These observations have raised many 
questions: what happens for flows without geometrical 
constraints such that fluctuations are of the same order 
of magnitude as the mean flow? Is the dynamo thresh- 
old strongly increased due to the lack of coherence of the 
driving flow 0, 0] or does the prediction made as if the 
mean flow were acting alone still give a reasonable order 
of magnitude [H? What is the nature of the dynamo bi- 
furcation in the presence of large velocity fluctuations? 
All these questions, and others motivated by geophysical 
or astrophysical dynamos 0], have led several teams to 
try to generate dynamos in flows with a high level of tur- 
bulence [13, EH ■ We present in this letter our first exper- 
imental observation of the generation of a magnetic field 
in a von Karman swirling flow of liquid sodium (VKS) 
for which velocity fluctuations and the mean flow have 
comparable kinetic energy and we discuss some of the 
above issues. 

The experimental set-up (see Fig. [IJ is similar to the 
previous VKS experiments [11| , but involves three modi- 
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FIG. 1: Sketch of the experimental set-up. The inner and 
outer cylinders are made of copper (in gray). The dimension 
are given in millimeter (left) and normalized by R (right). 
The 3D Hall probe is located either at point PI in the mid- 
plane or P2. In both cases, the probe is nearly flush with the 
inner shell. 



fications that will be described below. The flow is gener- 
ated by rotating two disks of radius 154.5 mm, 371 mm 
apart in a cylindrical vessel, 2R = 412 mm in inner di- 
ameter and 524 mm in length. The disks are fitted with 
8 curved blades of height h — 41.2 mm. These impellers 
are driven at a rotation frequency up to f2/27r = 26 Hz by 
300 kW available motor power. An oil circulation in the 
outer copper cylinder maintains a regulated temperature 
in the range 110-160 °C. The mean flow has the follow- 
ing characteristics: the fluid is ejected radially outward 
by the disks; this drives an axial flow toward the disks 
along their axis and a recirculation in the opposite direc- 
tion along the cylinder lateral boundary. In addition, in 
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the case of counter-rotating disks studied here, the pres- 
ence of a strong axial shear of azimuthal velocity in the 
mid-plane between the impellers generates a high level 
of turbulent fluctuations [H], [l3| ■ The kinetic Reynolds 
number is Re = KR 2 Vl/v, where v is the kinematic vis- 
cosity and K — 0.6 is a coefficient that measures the 
efficiency of the impellers [13]. Re can be increased up 
to 5 10 6 : the corresponding magnetic Reynolds number 
is, R m = K^ a aR 2 fl w 49 (at 120 °C), where ^ is the 
magnetic permeability of vacuum and a is the electrical 
conductivity of sodium. 

A first modification with respect to earlier VKS ex- 
periments consists of surrounding the flow by sodium at 
rest in another concentric cylindrical vessel, 578 mm in 
inner diameter. This has been shown to decrease the 
dynamo threshold in kinematic computations based on 
the mean flow velocity 14(. The total volume of liquid 



sodium is 150 1. A second geometrical modification con- 
sists of attaching an annulus of inner diameter 175 mm 
and thickness 5 mm along the inner cylinder in the mid- 
plane between the disks. Water experiments have shown 
that its effect on the mean flow is to make the shear layer 
sharper around the mid-plane. In addition, it reduces 
low frequency turbulent fluctuations, thus the large scale 
flow time-averages faster toward the mean flow. However, 
rms velocity fluctuations are almost unchanged (of order 
40 — 50%), thus the flow remains strongly turbulent [Hj]. 
It is expected that reducing the transverse motion of the 
shear layer decreases the dynamo threshold for the follow- 
ing reasons: (i) magnetic induction due to an externally 
applied field on a gallium flow strongly varies because of 
the large scale flow excursions away from the time aver- 
aged flow [l6| , (ii) the addition of large scale noise to the 
Taylor- Green mean flow increases its dynamo threshold 
, (iii) fluctuating motion of eddies increase the dynamo 
threshold of the Roberts flow [ItJ ■ 

The above configuration does not generate a magnetic 
field up to the maximum possible rotation frequency of 
the disks (H,/2tt = 26 Hz). We thus made a last modifica- 
tion and replaced disks made of stainless steel by similar 
iron disks. Using boundary conditions with a high per- 
meability in order to change the dynamo threshold has 
been already proposed [HI]. It has been also shown that 
in the case of a Ponomarenko or G. O. Roberts flows, 
the addition of an external wall of high permeability can 
decrease the dynamo threshold . Finally, recent kine- 
matic simulations of the VKS mean flow have shown that 
different ways of taking into account the sodium behind 
the disks lead to an increase of the dynamo threshold 
ranging from 12 % to 150 % f20|. We thought that using 
iron disks could screen magnetic effects in the bulk of 
the flow from the region behind the disks, although the 
actual behavior may be more complex. This last modi- 
fication generates a dynamo above R m ~ 30. The three 
components of the field B are measured with a 3D Hall 
probe, located either in the mid-plane or 109 mm away 



from it (PI or P2 in Fig. [T]). In both cases, the probe 
is nearly flush with the inner shell, thus B is measured 
at the boundary of the turbulent flow. Fig. [2] shows 
the time recording of the three components of B when 
R m is increased from 19 to 40. The largest component, 
By, is tangent to the cylinder at the measurement loca- 
tion. It increases from a mean value comparable to the 
Earth magnetic field to roughly 40 G. The mean values 
of the other components B x and B z also increase (not 
visible on the figure because of fluctuations). Both signs 
of the components have been observed in different runs, 
depending on the sign of the residual magnetization of 
the disks. All components display strong fluctuations as 
could be expected in flows with Reynolds numbers larger 
than 10 6 . 



20 



§10 




20 30 
time [s] 



FIG. 2: Time recording at PI of the components of the mag- 
netic field when the rotation frequency Q./2n is increased as 
displayed by the ramp below {Rm increases from 19 to 40). 

Fig. [3^ shows the mean values of the components (Bi) 
of the magnetic field and Fig. [3p their fluctuations Bi rms 
versus R m . The fluctuations are all in the same range (3 
G to 8 G, at 30% above threshold) although the corre- 
sponding mean values are very different. The time aver- 
age of the square of the total magnetic field, (B 2 ), is dis- 
played in the inset of Fig. [3^,. No hysteresis is observed. 
Linear fits of {B y ) or Bi rms displayed in Fig. [3] define a 
critical magnetic Reynolds number Rf n ~ 31 whereas the 
linear fit of (B 2 ) gives a larger value i?^ ~ 35. The latter 
is the one that should be considered in the case of a su- 
percritical pitchfork bifurcation. The rounding observed 
close to threshold could then be ascribed to the imper- 
fection due to the ambient magnetic field (Earth field, 
residual magnetization of the disks and other magnetic 
perturbations of the set-up). The actual behavior may 
be more complex because this bifurcation takes place on 
a strongly turbulent flow, a situation for which no rig- 
orous theory exists. The inset of Fig. [3b shows that 
the variance B 2 ms = ((B — (B)) 2 ) is not proportional 
to (B 2 ). Below the dynamo threshold, the effect of in- 
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duction due to the ambient magnetic field is observed. 
Brms/ ' (B 2 ) 1 / 2 first behaves linearly at low i? m , but then 
increases faster as R m becomes closer to the bifurcation 
threshold. We thus show that this seems to be a good 
quantity to look at as a precursor of a dynamo regime. 
In addition, we observe that it displays a discontinuity in 
slope in the vicinity of Rf n in an analogous way of some 
response functions at phase transitions or bifurcations in 
the presence of noise. Note however that the shape of 
the curves depends on the measurement point and they 
cannot be superimposed with a scaling factor as done for 
(B 2 ) versus R m in the inset of Fig. [3^. 

The above results are characteristic of bifurcations in 
the presence of noise. As shown in much simpler ex- 
periments, different choices of an order parameter (mean 
value of the amplitude of the unstable mode or its higher 
moments, its most probable value, etc) can lead to qual- 
itatively different bifurcation diagrams [2l[. This illus- 
trates the ambiguity in the definition of the order pa- 
rameter for bifurcations in the presence of fluctuations 
or noise. In the present experiment, fluctuations enter 
both multiplicatively, because of the turbulent velocity, 
and additively, due to the interaction of the velocity field 
with the ambient magnetic field. Finally, we note that 
both and are smaller than the thresholds com- 
puted with kinematic dynamo codes taking into account 
only the mean flow, that are in the range R c m = 43 to 150 
depending on different boundary conditions on the disks 
and on configurations of the flow behind them 14 , 2(| • 

The probability density functions (PDF) of the fluctu- 
ations of the three components of the induced magnetic 
field (not displayed) are roughly gaussian. The PDFs of 
fluctuations below threshold, i.e., due to the induction 
resulting from the ambient magnetic field, are similar to 
the ones observed in the self-generating regime. We do 
not observe any non gaussian behavior close to threshold 
which would result from an on-off intermittency mecha- 
nism [22j . Possible reasons are the low level of small fre- 
quency velocity fluctuations [23| or the imperfection of 
the bifurcation that results from the ambient magnetic 
field Q. 

Fig. |4]displays both the dimensional (see inset) and di- 
mensionless mean square field as a function of R m - (B 2 ) 
is made dimensionless using a high Reynolds number 
scaling @: (B 2 ) a p/(fx ((xR) 2 ) (R m - R c m )/R c m , where 
p is the fluid density. We observe that data obtained at 
different working temperature are well collapsed by this 
scaling (er decreases by roughly 15 % from 100 to 160 °C). 
The low Reynolds number or "weak field" scaling could 
also give a reasonable collapse of data obtained on this 
temperature range but the predicted order of magnitude 
for (B 2 ) would be 10 5 too small. 

Dissipated power by Ohmic losses is another important 
characterization of dynamo action. Our measurements 
show that, 30% above threshold, it leads to an excess 
power consumption of 15 — 20 % with respect to a flow 
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FIG. 3: a) Mean values of the three components of the mag- 
netic field recorded at PI versus R m (T = 120 °C): (A) -(B x ), 
(■) —{By), (•) (B z ). The inset shows the time average of the 
square of the total magnetic field as a function of R m , mea- 
sured at PI (•), or at P2 (*) after being divided by 1.8. b) 
Standard deviation of the fluctuations of each components of 
the magnetic field recorded at PI versus R m . The inset shows 
B rms /(B 2 } 1/2 . Measurements done at PI: (a) T = 120 °C, 
frequency increased up to 22 Hz; Measurements done at P2: 
(*) T = 120 °C, frequency decreased from 22 to 16.5 Hz, (□) 
T — 156 °C, frequency increased up to 22 Hz, (o) Q/2n = 16.5 
Hz, T varied from 154 to 116 °C, (0) ^/2tt = 22 Hz, T varied 
from 119 to 156 °C. The vertical line corresponds to R m — 32. 



driving power of the order of 100 kW. 

The effect of iron disks deserves additional discussion. 
A slight effect of magnetization of iron has been observed: 
the dynamo threshold during the first run was about 20 % 
larger than in the next runs for which all the measure- 
ments were then perfectly reproducible. However, no ef- 
fect of remanence that would lead to a hysteretic behavior 
close to the bifurcation threshold has been observed. De- 
magnetization of pure iron occurring for field amplitudes 
of the order of the Earth field, i.e., much smaller than 
the fields generated by the dynamo, the iron disks do not 
impose any permanent magnetization but mostly change 
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FIG. 4: The dimensionless quantity, {B 2 )no(oR) 2 / p is dis- 
played as a function of Rm for different working temperatures 
and frequencies (measurements done at P2 and identical sym- 
bols as in the inset of Fig. [3b) . The inset shows the same data 
in dimensional form B 2 versus rotation frequency for different 
temperatures. 



the boundary condition for the magnetic field generated 
in the bulk of the flow. This changes the dynamo thresh- 
old and the near critical behavior for amplitudes below 
the coercitive field of pure iron. It should be also empha- 
sized that the axisymmetry of the set-up cannot lead to 
Herzenberg-type dynamos [25| . In addition, these rotor 
dynamos display a sharp increase of the field at thresh- 
old and their saturation is mostly limited by the available 
motor power [25| . On the contrary, we observe a contin- 
uous bifurcation with a saturated magnetic field in good 
agreement with a scaling law derived for a fluid dynamo. 

The different mechanisms at work, effect of magnetic 
boundary conditions, effect of mean flow with respect to 
turbulent fluctuations, etc, will obviously motivate fur- 
ther studies of the VKS dynamo. A preliminary scan 
of the parameter space has shown that when the disks 
are rotated at different frequencies, other dynamical dy- 
namo regimes are observed including random inversions 
of the field polarity. Their detailed description together 
with experiments on the relative effect of the mean flow 
and the turbulent fluctuations on these dynamics are cur- 
rently in progress. 
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